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INTRODUCTION 

In these notes I used mainly the reviews / 1-3 written by Carsten Gundlach, 
I freely quoted in this paper sometimes even without quotation marks. If quota- 
tion makrks exist without references, these parts are from Gundlach reviews. This 
gentleman contributed to this field both with his original work and by the three 
reviews he wrote on this subject. 

I start by quoting Gundlach, from the first paragraph of still another paper 



that can be found in gr-qc/9604019 
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"In an astrophysical context, gravitational collapse normally starts from a star. 
This means that the initial data are almost stationary, and that they have a charac- 
teristic scale which is provided by the matter. Therefore, astrophysical black holes 
have a minimum mass, namely the Chandrasekhar mass. Abandoning the restric- 
tion to almost stationary data, or alternatively to realistic matter, one should be 
able to make arbitrarily small black holes. One may then ask what happens if one 
tries to make a black hole of infinitesimal mass by fine-tuning the initial data." 

We thus learn that for a stationary solution, a finite mass is necessary, but there 
may be surprises for the non-stationary solutions. Now I continue quoting Gundlach, 
but this time from the "Historical Introduction" part of his second review I 2 . 

"In 1987 Christodolou, who was studying the spherically symmetric Einstein- 
scalar model analytically, suggested to Matt Choptuik, who was investigating the 
same system numerically, the following question: Consider a generic smooth one- 
parameter family of asymptotically flat smooth initial data, such that for large 
values of the parameter p a black hole is formed, and no black hole is formed for 
small p. If one makes a bisection search for the critical value where a black hole 
is formed, does the black hole have a finite or infinitesimal mass?" "The idea was 
to start with a matter model that does not admit stable stationary solutions." 

"After developing advanced numerical methods for this purpose, Choptuik 
managed to give highly convincing numerical evidence that the mass is infinitesi- 
mal". Moreover he found two totally unexpected phenomena / 5 : 



The first is the now famous scaling relation 

M = c(p- P ,y i 

for the black hole mass M in the limit p close to but greater than p*. Choptuik 
found 7 approximately equal to 0.374 for all 1-parameter families of scalar field 
data. 

The second is the appearance of a highly complicated, scale-periodic solution for p 
close to . If we elaborate on this solution what is happening is the following. " For 
marginal data, both supercritical and sub-critical, the time evolution approaches 
a certain universal solution which is the same for all the one-parameter families 
of data. This solution is an "intermediate attractor" in the sense that the time 
evolution first converges onto it, but then diverges from it eventually, to either 
form a black hole or to disperse. This universal solution, also called the "critical 
solution", has a curious symmetry." The critical solution 0*(r, t) is the same when 
we rescale space and time by a factor e A : 

4>* = <P*(e A r, e A t). 2 

The logarithmic scale period of this solution, A approximately equal to 3.44 is a 
second dimensionless number coming out of the blue. 

As a third remarkable phenomenon, both the " critical exponent" and " critical 
solution" are "universal", that is the same for all one-parameter families ever in- 
vestigated independent of the shape of the initial data, whether it is a gaussian, a 
hyperbolic tangent or the two other forms studied in reference 5. For a finite time in 
a finite region of space, the space-time generated by all near-critical data approaches 
one and the same solution. This universal phase ends when the evolution decides 
between the black hole formation and dispersion. The universal critical solution is 
approached by any initial data that are sufficiently close to the black hole threshold, 
on either side, and from any 1-parameter family. 

Similar phenomenon to Choptuik's results were quickly found in other systems 
too, even in gravitational waves ( which can form black holes in the absence of 
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matter) suggesting that they were limited neither to scalar field matter nor to 
spherical symmetry. 

The cylindyrically symmetric case was studied by A.M. Abrahams and C.R. 
Evans / 6 with similar results, having A approximately equal 0.6 and 7 approxi- 
mately equal to 0.36. For a perfect fluid with equation of state p = | I 1 , the 
"critical solution" is not discretely, but continuously self-similar (CSS) in spherical 
symmetry, since the universal attractor has a different symmetry. 7 is found to be 
0.36 once more. Other models were also studied with different values for 7. For 
models in which fermions are coupled to gravity, similar phenomenon is found / 8 . 
From these examples it was clear that 7 was not universal with respect to different 
kinds of matter, but only work with respect to the initial data for any one matter 
model. 

When similar systems were studied in lower dimensions, namely in 2 + 1, an- 
alytical results were also possible. When the system was studied numerically by 
Choptuik and Pretorius / 9_13 ? one obtained what seemed to be irrational indices. 
When Birmingham / 14 studied the "properties of a massive scalar field in the back- 
ground geometry of the 2+1 dimensional black hole" , he could solve the wave equa- 
tion exactly in terms of hypergeometric functions . The remarkable feature was 
that he could evaluate one of the Choptuik scaling parameters as the sum and dif- 
ference of the outer and the inner horizons and show that it was given by a rational 
number, namely \ . This reminds us that the 2+1 dimensional black hole, so called 
BTZ black hole, exists in a space which is 'related' to the anti-de Sitter space. The 
latter is related to a conformal model / 15_17 in one less dimension by the Malda- 
cena correspondence Z 18 . The emergence of hypergeometric functions and rational 
exponents are signatures of "conformal invariance"/ 19 ' 20 . If we consider only the 
relation of hypergeometric functions to conformal symmetry also found in / 21_23 . 
We will comment on lower dimensional solutions later. 

If we go back to 3 + 1 dimensions , we see that interesting phenomenae occur 
in self-similar solutions. General related phenomena divide into two categories, the 
continuosly self-similar solutions like those found by Evans and Coleman/ 7 where 
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the metric coefficients are of the form 



9i*v(t, x l ) = e T g* fl „(x l ), 3 

and discretely self -similar solutions, found by Choptuik/ 5 to where the metric 
coefficients are of the form 

g^TiX*) =e- 2T g* ll „(T,x l ). 4 

The right hand side of this equation is given as 

g*^ (r,x l ) = g * flu (t + A,x l ). 5 

Here r is the negative of the logarithm of a space-time scale, r = —log{r 2 + t 2 ), and 
Xi are dimensionless. 

Work in related fields still continue. A last paper I found by Choptuik et al is 
given in |gr-qc/0305003| Z 24 . We must first study the first work on this subject. 



THE SPHERICALLY SYMMETRIC CASE WITH SCALAR FIELD 

Choptuik / 5 studied a massless scalar field (p minimally coupled to the gravi- 
tational metric. 

= 8ttT^ 6 

where LHS side has the Einstein tensor and the RHS has the stress-energy tensor 
of the scalar field. The scalar field obeys the Klein-Gordon equation written in this 
background. 

He chose the metric as 



ds 2 = -a 2 (r, t)dt 2 + a 2 (r, t)dr 2 + r 2 rf0 2 , 



where dfl 2 = d0 2 + sin 2 9d(p 2 is the metric on the unit 2- sphere. He chose t 
orthogonal to r to have no drdt cross term in the line element. To fix the coordinate 
completely, he chose a = 1 at r = 0, so that t is the proper time of the central 
observer at r = 0. 



He used auxiliary variables 



n = -4>, t 8b 

a 



to make the wave equation a first order system: 



$ , = ( -n 



t 

a J ,r 



r z V a ) , r 
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In spherical symmetry there are four algebraically independent components of 
the Einstein equations. One of them can be written in terms of derivatives of the 
others. The others read 

^ + - 2yrr(n 2 + $ 2 ) = 0, 11 

a 2r 

a r a r a 2 — 1 „ 

= 0, 12 

a a r 

-± - 4yrr$n = 0. 13 

a 

Choptuik chose to use the first two equations with no time derivatives on the metric 
parameters. The scalar field evolved forward in time, while a and a were calculated 
from this result at each new time step, by explicit integration over r starting at 
r = 0. The third equation obeyed automatically. He took a = 1 at r = to make 
the spacetime regular at that point. These equations do not have an intrinsic scale, 
and scale transformation on t and r transform one solution into another. 

There is no gravitational radiation for the spherically symmetric case. The 
degrees of freedom of the gravitational field depends on the matter fields. Thus the 
free data for the system will be given in terms of two functions n(0, r) and $(0, r); 
the presence of gravity does not change this fact. 

Choptuik studied 1-parameter families of such data by evolving the data for 
many values of the parameter p, which he chose to be the amplitude of the Gaussian 
he took for n(r, 0). He could have taken the width or the center of the Gaussian as 
well with the same result. He took $(0, r) = 0. 
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Christodolou / 25 > 26 has shown that for spherically symmetric systems, data 
sufficiently weak in a well-defined way evolve into a Minkowski-like spacetime, and 
that strong data forms a black hole Z 27 . What Choptuik was looking at what 
happens in between, where the conditions of neither theorem apply. 

He found that he could make arbitrarily small black holes by fine-tuning the 
parameter p ever closer to the black hole threshold. One point was important. 
Nothing singled out the black hole threshold in the initial data. One could not tell 
whether one set of data would form a black hole whereas another set infinitesimally 
close to this one would not. One only had to let the scalar field evolve in time 
for a sufficiently long time. At the end he showed that one could form black holes 
as small as 10~ 6 times the ADM mass of the spacetime. The power law scaling 
was obeyed over six orders of magnitude of masses. He found that the scalar wave 
contracts until curvature values of order 10 12 are reached in a spacetime region of 
size 10 -6 before it starts to disperse. He found that all near-critical spacetimes, 
for all families of initial data, looked the same in the intermediate region. They 
approximate one universal spacetime, called the critical solution 

Z*(kr,k(t-U)). 14 

The accumulation point and the factor k depend on the family, but Z* does not. 
This universal solution has the property that 

Z*(r,t) = Z*(e nA r,e nA t) 15 

for all integer n and for A = 3.44... Here Z stands for the metric coefficients and 
the scalar field. 

All critical points that have been found in black hole thresholds so far have 
an additional spacetime symmetry . They are either time- independent or scale 
invariant. These two kinds exhibit different phenomena, the former giving rise to 
continuous and the latter to discrete self similar solutions. 

Those that are time independent have a mass gap. The critical solution is by 
definition an attractor of co-dimension one, which means that it has precisely one 
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unstable perturbation, with Ao greater than zero. All other perturbations decay, 
having the real part of their eigenvalues Aj less than zero. Intuitively, these critical 
solutions can be thought of as metastable stars , and they typically occur when the 
field equations set a mass scale in the problem. 

The second type critical solutions occur when there is no scale in the field 
equations, or when this scale is not dynamically relevant. Many systems, such as 
a massive scalar field, show both types of critical phenomena in different regions of 
the space of initial data. 

For the scale-invariant system too, this symmetry comes in a continuous or 
a discrete version. A continuously self-similar (CSS) solution is invariant under 
an infinitesimal rescaling of both space and time. A discretely self-similar (DSS) 
solution is invariant only under rescaling by a particular finite factor, or its integer 
powers. They are independent of a suitable scale coordinate. They give rise to 
power-law scaling of the black-hole mass at the threshold. 

Similar phenomena occur in statistical mechanics where self-similar solutions 
near a second order phase transition are important. At the critical point between a 
gas and a liquid, and otherwise clear liquid becomes opaque, due to density fluctu- 
ations appearing on all scales up to scales much larger than the underlying atomic 
scale, and including the wavelength of light. The density discontinuity vanishes as 
a non-integer power: 

Pliquid ~ Pgas = k(T* — T) 7 . 16 

Similar behaviour is seen in a ferromagnet where the spontaneous magnetization 

goes as 

|m| = c(T* — T) 7 . 17 

In both of these expressions T* is the critical temperature where the respective 
"order parameter" vanishes. At this point all the quantities that are dimensional 
in the model are 'irrelevant', and the system is approximately scale-invariant. 

Further work was done on the interpretation of the initial Choptuik solution. 
A most recent study of the global structure of this solution was done by Martin- 
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Garcia and Gundlach / 28 where "this solution, formed in the gravitational collapse 
of a scalar field resulting in a naked singularity, was found to be regular at the center 
to the past of the singularity, and regular at the past lightcone of the singularity. At 
the future lightcone of the singularity, which is also a Cauchy horizon, the curvature 
is finite and continuous but not differentiable. To the future of the Cauchy horizon 
the solution is not unique, but depends on a free function ( the null data coming 
out of the naked singularity). There is a unique continuation with a regular center 
(which is self-similar). All other self-similar continuations have a central timelike 
singularity with a negative mass Z 28 ." 

SOLUTIONS IN 2+1 DIMENSIONS 

In 2+1 dimensions, where we do not have to stick only to numerical calculations, 
analytical solutions may be easier to obtain and tell us more about the underlying 
mechanism of this phenomenon . General relativity in 2 + 1 dimensions, however, 
is quite different from that in 3 + 1 dimensions, because the degree of freedom 
of the graviton is given by d — 3 in any number of dimensions. This count does 
not allow any local degrees of freedom to the graviton in 2 + 1 dimensions. One 
could study models only with interesting global properties if we study pure gravity. 
Furthermore, in 2 + 1 dimensions the Riemann tensor does not have any other 
components as those given by the Ricci tensor, i.e. the Weyl tensor which gives 
half of the components of the Riemann tensor in 3 + 1 dimensions is identically 
zero. Another tensor, called the Cotton tensor, which exists in 2 + 1 dimensions 
is also identically zero when we use the Einstein-Hilbert action for gravity. This 
fact makes it impossible to have a Ricci flat space with non-trivial Riemann tensor 
components. Due to this fact, there are no gravitational waves. Blackholes can 
exist only in anti de Sitter spaces Z 29 , rather than asymptotically flat spacetimes. 
This solution is called the BTZ solution. 

Quoting from Gundlach I 1 "The negative cosmological constant of the AdS 
space locally introduces a mass scale into the field equations. It also changes the 
the global structure of spacetime. Null infinity becomes a timelike surface. A null 
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ray can go out from the center to null infinity and return while a finite proper time 
passes at the center. The only consistent boundary conditions for the massless wave 
equations are the Dirichlet boundary conditions. This means that all outgoing scalar 
waves are reflected back to the center in a finite time." In a later article / 9 Choptuik 
et al. studied " Gravitational collapse in 2 + 1 dimensional AdS spacetime" . They 
find what looks like a CSS solution, of the second type ( with mass vanishing at the 
critical point). "In this respect, one should note one thing. Since the cosmological 
constant introduces reflecting (Dirichlet) boundary conditions for the scalar field, 
all of this matter must eventually fall into the black hole. In this sense, the black 
hole mass is simply the asymptotic mass of the spacetime, and there is no black hole 
threshold. In numerical calculations, however, one uses initial data with compact 
support. One can define the black hole threshold as the formation of an apparent 
horizon before the scalar waves are reflected for the first time. An apparent horizon 
is the surface where the gradient of r is null. If the cosmological constant is negligible 
and the mass M is much less than unity, the dynamics is approximately scale 
invariant and second type solutions may be found." I 1 

Pretorius and Choptuik / 9 found M going as (p — p*) 27 where 7 = 1.20 ap- 
proximately. Husain and Olivier / 10 studied the similar system and found 7 = .81. 
Garfinkle Z 13 studied this case analytically by setting the cosmological constant to 
zero, and found an index close to that of Pretorius and Choptuik. Garfinkle and 
Gundlach / 30 calculated the perturbation spectrum of Garfinkle's solutions in closed 
form with the assumption that the perturbations are analytic. They found instead 
of one growing mode, three growing modes. The dominant one gives the index for 
an unstable mode to be equal to | which is in good agreement with the Pretorius- 
Choptuik work. All this work has uncertainties , though. Another index found in 
reference 30 is |. Hirschmann et al / 31 use a constraint to suppress some of these 
modes and find the index 7 equal to four. 

Birmingham and Sen / 32 considered the formation of a black hole from the 
collision of two point particles of equal mass in 2 + 1 dimensions with a negative 
cosmological constant, which would have the metric given below in the absence of 
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the scalar particles. 



ds 2 = - (-M + -^—-] dt 2 + (-M + -^—-] dr 2 + r 2 d9 2 . 18 

V (-a); V (-a); 

Here M takes values from —1 to infinity. M = — 1 is the anti-de Sitter space. 
— 1 < M < results in point particle naked singularities at the center. When M 
equals or exceeds zero, we have black holes. 

Birmingham and Sen found the index to be \ for the case when the phase space 
is two dimensional. Peleg and Steif Z 33 studied the collapse of a dust ring where 
the space of initial data is two dimensional and found the same value for the index. 
In the Birmingham / 34 construction the presence of two Virasoro algebras, so the 
relation to conformal symmetry in two dimensions is evident. He goes further and 
expands a massive scalar field in the background metric of a rotating black hole: 

,2 ( „ r 2 J 2 \ ;2 / „ r 2 J 2 \ _1 l2 <i( m Jdt ^ 2 
ds 2 = - —M + — — + — r dt 2 + -M + — — + — dr 2 +r 2 d6 - 



(-A) 4r 2 J V (~ A ) 4r2 / V 2r2 

19 

In this background the scalar wave equation can be solved by an ansatz of the form 

$ = R{r)e- iu}t e im(t> . 20 
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After a change of variables to z = ^ 2 _^j" , where r+ and r_ denote the outer and 
inner horizon radii, the points when (— M + -^-^ + = 0, the radial equation 
is of the hypergeometric type, with solution R(z) = z a (l — ^f-Fi^fr, c, z). The 
implementation of the vanishing boundary condition at infinity, fixes the constants 
a, 6, c, a, [3 in terms of the parameters of the metric, which can be used to find the 
index for the mass, after a detour . 

As stated above, the presence of hypergeometric functions among the solutions 
reminds us of the other places where they occur, always related to the presence 
of conformal symmetry in the problem. We should note that we get rational in- 
dices when we perturb around the exact solution which is related to a model with 
conformal invariance 15_17 . 

In still another paper Birmingham et al. Z 35 show how the 2 + 1 black hole 
gives rise to: 
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a) no hair theorem: The Euclidean BTZ black hole has to topology of a torus which 
is parametrized only by two real parameters, which makes one conclude that at 
most two parameters, mass and angular momentum are sufficient to describe the 
system. At this point also note how one can not construct a charged rotating black 
hole directly, but have to use an indirect method Z 36 . 

b) holography: BTZ black hole is a holographic manifold, since the three dimen- 
sional hyperbolic structure is in 1-1 correspondence with the parameters of the 
two-dimensional genus one boundary. 

c) entropy: The Bekenstein-Hawking entropy is determined by the Teichmiiller space 
on the boundary. 

d) Maldacena Conjecture Z 18 : The three dimensional hyperbolic space is related to 
a SL{2, Z) family of solutions. 

We will continue this review with a summary of some work / 37 in 2 + 1 di- 
mensions where some of the above mentioned phenomena is displayed a little more 
clearly. We will show, how in this overly simplified case, the naked singularity is 
avoided. 

TOY MODEL 

We start with the equations given by Pretorius and Choptuik / 9 for a metric 
given by 

ds 2 = ^—(dr 2 - dt 2 ) + tan 2 (r)e 2B ^d9 2 . 21 
cos z (r) 

Here the full space is mapped into the interval < r < ^. The cosmological constant 
is chosen equal to —1 and r is scaled so that it is a dimensionless parameter. We 
take the static and the spherically symmetric case, where A, B are functions of r 
only. Since we are treating the spherically symmetrical case, 9 independence is 
justified. Here, as in all static cases, we expect to have a finite mass. We are using 
this example only to show how the singularity is formed, and how it is avoided. 

If a scalar, static particle is coupled to the static metric, the Einstein equations, 
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with the cosmological constant read / 9 

(1-e 



2A\ 



A,rr + 1 TTT~ + 27r ^ r = 0, 22 

cos A [r) 

4 \ 2(1 -e 2A ) 



B rr + B r B r + , , r H 2 , >. = 0, 23 

\ sm(2r) J cos z (r) 

x, (l + cos 2 (r))\ 2A r (l-e 2A ) n j2 rt rtJ 

\ szn(2r) / sin(2r) cos^(r) 

[tanre B '4> )r \,r = 0. 25 
This system of equations have the set of solutions 

0=0, ,4= -log(sinr),B = log( =-). 26 

2sin 2 r 

The expressions for A and (p were given by Pretorius and Choptuik . We see that 
our original domain, < r < ^ is halved. Due to the singularity in £?, we can use 
only the region where < r < j. 

We perform a simple perturbation expansion at this point treating the above 
set of solutions as the zeroth order term. 

A = A + eAt + e 2 A 2 + , 27 

B = B + eB 1 + e 2 B 2 + , 28 

= e0i + e 2 2 + 29 

Equation (25) gives us 4>i,r = Ctan2r where C is a constant of integration. We 
choose C = \. The equation ( 22 )reduces at this order to 

A 8A -0 30 

which can be reduced to an equation of the hypergeometric type. A simple calcu- 
lation shows that we have a special form of the hypergeometric equation, yielding 

sin- 1 (2r) 2 F 1 (-l/2, -1/2| - l/2|sm 2 (2r)) 31 
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for one solution, and 

sm 2 (2r) 2 F 1 (l,l|5/2|sm 2 (2r) 32 
for the other. These special forms of the hypergeometric function reduce to 

A\ = cot(2r), 33 

A\ = 3(1 - 2rcot2r) 34 
respectively. For the former solution A\ we get two solutions 

B\ = cot2r + tan2r, 35 

B 2 = -6rBl 36 

By extending the perturbative analysis, for A we get a differential equation, again 
of hypergeometric type , with an inhomogenous term, for A which can be easily 
integrated. We are rather interested in the solution for 0, to study its singularity 
structure. To second order in e 

12r — 3sm(4r) 
02 ' r = 2««2(2r) ' 

In the expression found for 0, there is no sign of singularity at r = 0. 

Our results may be interpreted better if we transform our original coordinates 

to 

f = tan(r)e B{r) . 38 
Then our metric is transformed into 

_9 dv 2 

ds 2 = —i—M + f 2 )dt + + rW, 39 

(-M + r 2 ) 

where 

—M = e 2{B ~ A) [sec 2 (r) + 2tan(r)B, r + sin 2 (r)B 2 r ] - tan 2 (r)e 2B 40 

as given in the Pretorius-Choptuik reference . Note that here M is r dependent. 
In the original case / 9 , it is also time dependent. Using our zeroth order solution 
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for A, we find f — cot(2r) + 0(e). This transformation maps our original domain, 
< r < f into oo > r > 0. 

For the zeroth order solution — M = 1, a constant. We see that our solution 
corresponds to the AdS solution, known for this system. If we use the solution 
B = log( 2 ^°J^ r z r ), where < a < 1,, we get the solutions corresponding to conical 
singularities Z 38 . 

The BTZ solution is obtained if we set a = — i which necessitates a reparametriza-| 
tion using hyperbolic functions. In the parametrization of the metric , eq.(21), hy- 
perbolic cosine and hyperbolic tangent functions replace trigonometric cosine and 
trigonometric tangent respectively. The presence of i in the expression in the ex- 
pression for B introduces the necessary minus sign for — M, eq.(39), while replacing 
tan(r) by tanh(r) in the metric, eq(21), retains the original signature. 

We, then, use our perturbative solutions in the presence of the scalar field. At 
first order in e, we find we have two solutions. If we take the solution set A\, B\, 
we find 

-M = 1 + 32e[= -r 3 ]. 41 
r 

This expression diverges both at f = and r = n/4 in an undesired fashion, so is 
discarded. For the latter solution set, A\,B\, 

-M = 1 -96e[(l + r 2 ) + cot' 1 (r)(-r 3 + -)], 42 

an expression which diverges at the origin, but this time with the correct sign. 
At r = 0(r = 7r/4), — M diverges to minus infinity. We excise the domain when 
— M + r 2 equals zero. 

When the scalar field is calculated in terms of the new variable, we find 

cj) = 1/2/0(7 ( ( ] J • 43 

For r = 0, the scalar field cf> is proportional to log(r), as pointed out by Garfinkle 
/ 13 and Burko Z 39 . When r goes to infinity, r = 0, goes to zero. The similar 
behaviour persists at second order in e, where 

02,r = -3=-cor 1 (r)(^ = ti). 44 
r r 
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We see that the divergence is severer when r goes to zero. 

The fact that — M goes to minus values as r approaches zero, signals the pres- 
ence of a black hole around the origin. We can not tolerate — M + f 2 being null. 
We excise the space at the value of f where — M + f 2 equals zero. There is only one 
root of the equation 

—M + r 2 = 0, 45 

giving the approximate condition 

r > 24ne (l - (24) Ve 2 )) +0(e 4 ). 46 

This will also prevent the curvature singularity which will occur at r = 0. 

We have the scalar curvature made out of two parts , the finite part corre- 
sponding to the AdS solution and a singular part coming from the perturbative 
solution. 

^=-6 + ^— 2 47 
r 1 + r 

where R is the curvature scalar given by Pretorius and Choptuik. 

To detect whether black-hole is formed or not, we use a second test given by 
Pretorius and Choptuik, by checking the condition for "trapped surfaces". In this 
reference, " trapped surface" is defined " to be surfaces where the expansion of the 
outgoing null curves normal to the surface is negative " . The condition for this to 
happen is given at the same place as 

S = 1 + sin(r)cos(r)B jr < 0, 48 

as applied to our case. We study whether this constraint is satisfied for our solutions. 
To zeroth order in e this condition reads 

(1 + (r 2 ) 1 / 2 

2S = - { 1 } . 49 

r 

We see that "the surface is trapped" for r = (r = 7r/4). For r = i.e. r 
approaching infinity, we get S = — which goes to zero from below. If we use the 
first order solution, however, the correction changes the situation in the first case, 
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whereas it does not change the result in the second case.. For the first solution, we 
have 

S = 1 + sin(r)cos(r)B tr = 

For r/4 < e there is no "trapped surface. We can come all the way to the curvature 
singularity. This solution is discarded also for the above reason. If we use the second 
solution, however, we get 

S = -(l + (r)Y /2 (h 
r 



-3e 



(1 + (r) 2 )" 1 / 2 ((+^+ =) + cat-\r) {-r 2 + ^ 



51 



As r goes to zero, the surface is trapped. The particle is not allowed to come close 
to the coordinate independent curvature singularity at the origin. 

After we find the expressions given above for the scalar field and the metric, we 
use these approximate solutions, eq.s (35,36,43,44) and check whether they satisfy 
the set of exact equations, eq.s(22-25).The figures 1 to 4 of this reference Z 37 show 
that although our approximate solutions do not satisfy the equations as r goes to 
zero, the behaviour as r goes to infinity of these figures clearly show that asymptot- 
ically these solutions tend to be exact. We thus see that we can obtain approximate 
solutions to the equations of motion and the constraint equations, which approach 
exact solutions in the asymptotic region. Our tests, mainly the fact that — M + f 2 
goes through zero tells us that a black hole is formed as r goes to the origin (figure 
5 of the reference). We excise the space at the point where —M + f 2 equals to zero. 
At this point note that our first set of approximate solutions, eq.s (33,35) do not 
give the correct asymptotic behaviour, so are discarded. 

From these figures we see that our approximate solution is no longer reliable 
as f goes to zero. Since we excise our space in this region and the solution has the 
correct asymptotics as r goes to infinity, we think that the message our approximate 
solution conveys, i.e. the presence of a black hole at the origin, is correct. 
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Final Remarks 

Here I tried to give a taste of a subject Matt Choptuik started a decade ago. 
there is a large literature built around this work nowadays. The original Choptuik 
paper 5 has over three hundred citations in ten years, which may not be a lot 
compared to a paper on a particle physics experiment, but is considerable, when 
compared to the citations of one of the best known reviews in general relativity, 
namely the Eguchi, Gilkey, Hanson paper / 40 (670 citations in twenty three years). 
This topic is considered one of the areas where most of the work numerical relativity 
is done. 
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